Geometric (Berry) phases and statistical transmutation in the two-dimensional 

systems of strongly correlated electrons 
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Focusing on the hole-doped two-dimensional systems of strongly correlated electrons, we examine 
geometric phases acquired by electronic wave functions as a result of polaron transport around a 
closed loop. For this study we apply the Lanczos exact diagonalization method to Holstein-Hubbard, 
Holstein-iJ, and Holstein-iJ^ models in order to reveal various aspects of geometric phases. We 
demonstrate that transverse spin fluctuations are responsible for the generation of nontrivial geo- 
metric phases. From the exchange symmetry of polarons we flnd that the statistical transmutation 
of polarons occurs depending on the strength of electron correlations. 
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I. INTRODUCTION 

Since the discovery of high temperature superconduc- 
tivity in copper oxides, attempts have been made to 
identify the nature of charge carriers in these materi- 
als. There are some experimental evidences that pola- 
ronicj-|Charge carriers exist in rbpith the superconducting 
stateEl and the normal statefl'cl Photoinduccd absorp- 
tion measurements in La2Cu04 and YBa2Cu307_5 indi- 
cate that self-localized structural distortions are present 
around photoexcited charge carriersJj A recent experi- 
ment demonstrated that Sr-induced local lattice distor- 
tions occurs in La2-xSrxCu04 in association with holes 
donated by the Sr atoms.H In order to investigate the na- 
ture of the polaronic charge carriers, we study geometric 
phases acquired by the electronic wave functions after the 
transport of a polaron around a closed loop. 

A quantum object (a subsystem of fast degrees of free- 
dom) acquires a geometric (Berry) phase in an environ- 
ment where a parameter (a subsystem of slow degrees 
of freedom) is slowly transported around a closed cir- 
cuit in a parameter space.B Thus quantum subsystems 
with two widely separated energy scales (between sIqi 
and fast degrees of freedom) exhibit geometric phases B 
In this paper we investigate the geometric phases ac- 
quired by the electronic wave functions as a result of 
polaron transport around a closed loop in a hole-doped 
two-dimensional system of strongly correlated electrons, 
e.g., the copper oxide plane in high Tc superconductors. 
In addition, we investigate the statistical transmutation 
of polarons by studying the exchange symmetry of po- 
larons varying with the strength of electron correlations. 



II. VARIOUS GEOMETRIC (BERRY) PHASES 
BASED ON MODEL HAMILTONIANS 

We first study the geometric phases varying with the 
strength of antiferromagnetic electron correlations and 
of electron-phonon coupling. For this study we introduce 
the Holstein-Hubbard model Hamiltonian for the hole- 



doped two-dimensional systems of square lattice,Llll3 
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Here cj^ (c^o-) is the creation (annihilation) operator of 
an electron with spin cr at a copper site i (i = {ix,iy)), 
and n^^ — cl^c.^^ is the electron number operator, xi and 
yi denote the displacements of oxygen atoms in the unit 
cell of Cu02 plane.EilO t represents the electron hopping 
integral; U, the on-site Coulomb repulsion energy (cor- 
relation strength); g, the electron-lattice coupling con- 
stant and K, the spring constant. In our calculations 
the magnitude of the Holstein coordinate Uj is normal- 
ized to unity, \u^\ = 1 when a polaron is formed around 
copper site i as a result of hole doping. The energy of 
the polaron is then given by —g. 

Now for the study of geometric phases varying with 
the Heisenberg coupling constant J and the electron- 
phonon coupling constant g, Hamiltonians of interest are 
the Holstein-iJ model,BE3 
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and the Holstcin-t Jz model. 
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Here 



= Cj^(l — "i _ct) is the electron annihilation 
operator at site «, which excludes double occupation. 



^i — l^ia'^ aR^iis ^^ ^^^ electron spin operator. The 
Holstein-t Jj Hamiltonian takes into account only the spin 
^-component interaction. 

Thfipiattice distortions are treated in the adiabatic 
limit .Q'tHl The electronic wave functions (the fast degree 
of freedom) varying with lattice distortions are obtained 
from the, application of Lanczos exact diagonalization 
methodt^ to a tilted— j/IO x \/T0 lattice with periodic 
boundary conditions .113 The propagation of the local lat- 
tice distortion (polaron hopping) is described by 
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with r being the dimensionless time lapse with < r < 1 
for polaron hopping between adjacent lattice sites. For 
instance, the lattice distortion occurs at copper site A at 
r = and moves to an adjacent copper site i? at r = 1. 
In this manner, a closed path in the parameter space {u^} 
can be defined. We compute the geometric phase factors. 
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acquired by the electronic wave function after po- 



laron transport around a closed path C during time T, 
by usingtil 
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Here |n(u(rfc))) is the electronic eigenfunction for a given 
lattice distortion u(rfe) at time Tfe. Tfe is a fcth discretized 
time lapse between the initial time Tq = and the fi- 
nal time (period) r„ = T. Eq. (^ is derived under the 
condition that \n{u{Tk))) is a single- valued-complex wave 
function, that is, \n{u{0))) = \n{u{T)))Bk3 On the other 
hand, all the electronic wave functions obtained from 
the diagonalization of the model Hamiltonians above are 
real- valued. In such cases, geometric phases a«pjea,r in 
the form of double- valued real wave functions .Il3"ll3 Us- 
ing the local gauge transformation in parameter space 
w, the double-valued real wave function, \n{u))r can be 
converted into a single-valued complex wave function 
\n{u))c, 
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with 0{u), a differentiable function and 0{u{T)) — 
9{u{Q)) = TT (see Appendix for details). By introducing 
the above expression into Eq. (H), one can now obtain 
the geometric phases of interest. 



III. COMPUTED BERRY PHASES AND 
STATISTICAL TRANSMUTATION 

First, by using the Holstein-Hubbard model we con- 
sider the Berry phase acquired by the electronic wave 
function as a result of polaron hopping around a closed 
path. Two types of closed paths, a triangular path and a 
square path are displayed in Figs. |l|(a)-(b). For simplic- 
ity, for the next nearest neighbor hopping we chose the 



identical value of hopping integral to the one used for the 
nearest neighbor hopping. The electronic wave function 
is predicted to gain a geometric phase by -k (the geo- 
metric phase factor of —1) for polaron hopping around 
the smallest possible closed path, that is, the triangu- 
lar path [Fig. |l|(a)]. On the other hand, the electronic 
wave function with the polaron transported around the 
square path [Fig. |^(b)] is predicted to gain a phase an- 
gle by 27r, thus leaving the phase factor of the electronic 
wave function unchanged. The geometric phase of 27r 
can be understood from the decomposition of the square 
path into two triangular paths, as shown in Fig. S; the 
geometric phase factor of -1-1 for the square path results 
from the product of two identical geometric phase factors 
of —1 for the decomposed triangular paths. 
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FIG. 1. Geometric phase factors based on the Hol- 

stein-Hubbard model for (a) a triangular loop and (b) a square 
loop. The horizontal axis denotes the Coulomb correlation 
strength U and the vertical axis, the electron-lattice coupling 
constant g. The energy unit is i = 1. 
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FIG. 2. The decomposition of a square path into two tri- 
angular paths. 

By using both the Holstein-i J model and the Holstein- 
tJz model, we now examine the geometric phases which 
result from the transport of a single polaron. They are 
shown in Figs. ^(a)~(b) and Figs, ^(a)-(b). The pre- 
dicted geometric phase factor for the case of the tri- 
angular path [pig. 0(a)] is in agreement with that of 
Schiittler et a/.ll3 For large U values, say, U — 8t both 
the Holstein-t J model and the Holstein-Hubbard model 
yielded identical geometric phase factors for both the 
triangular path and the square path, as are shown in 
Figs. |l](a)-(b) and Figs. Q(a)-(b). Such direct compari- 
son between the two model Hamiltonians is valid only for 
large U values. This is because the t-J model is equiva- 
lent to the Hubbard model Hamiltonian corresponding 
to large U values. On the other hand, the Holstein- 
tJz model yields a trivial geometric phase factor of -1-1 
for the triangular path [Fig. 0(a)], while the Holstein- 
tJ model Hamiltonian gives rise to the nontrivial geo- 



metric phase factor of —1. It is of note that only the 
t-J (but not the t-Jz) Hamiltonian allows transverse 
spin fluctuations. Thus we discover from comparison of 
the two Hamiltonians that the generation of the non- 
trivial geometric phase factor of —1 is caused by the 
transverse spin fluctuations (or spin flip-flop fluctuations) 
which comes from interactions involving the x- and y- 
component (but not the z-componcnt) spins, that is, 

■^ 22{ij) WixSjx + SiySjy) = 2 "^ 2_/(y) [Si+Sj- + Si-Sj+). 

The nontrivial geometric phases are also predicted by 
the Holstein-Hubbard model Hamiltonian, in agreement 
with the Holstein-iJ model Hamiltonian, as is shown in 
Fig.g(a). 



FIG. 3. Geometric phase factors based on the Holstein-t J 
model for (a) a triangular loop and (b) a square loop. The 
horizontal axis denotes the Heisenberg exchange interaction 
J and the vertical axis, the electron-lattice coupling constant 
g. The energy unit is i = 1. 



(f/^6). It is thus proper to state that polarons behave 
as hard core bosons in the case of weakly and strongly 
correlated electrons, while they act as fermions in the 
intermediate region of electron correlations. Such statis- 
tical transmutation of polarona_can be mapped into the 
Jordan- Wigner transformationEj for the two-dimensional 
lattice, 



iS t 



(7) 



We find that the predicted phase angle S depends on U, 
the strength of electron correlation (Coulomb repulsion). 
That is, (5 = TT for small and large U values, indicating 
that polarons behave as bosons and 6 = for intermedi- 
ate U values, indicating that they act as fermions. 
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FIG. 5. Geometric phase factors based on the Hol- 

stein-Hubbard model for two-polaron exchange around a 
square path. The numeral label for the path indicates the 
consecutive order of polaron hopping. The axes are the same 
as in Fig. hi. 
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FIG. 4. Geometric phase factors based on the Holstein-t J^ 
model for (a) a triangular loop and (b) a square loop. The 
horizontal axis denotes the Ising interaction Jz and the verti- 
cal axis, the electron-lattice coupling constant g. The energy 
unit is f = 1. 

Now we investigate the exchange symmetry of two po- 
larons based on the Holstein-Hubbard model. The nu- 
meral label for the square path in Fig. g indicates the con- 
secutive order of polaron hopping. Interestingly enough, 
there exist both trivial and nontrivial geometric phase 
factors. After the two polarons are exchanged, the elec- 
tronic wave function is seen to acquire a trivial geometric 
phase factor of -f 1 for the case of weak electron correla- 
tions (C/ $ 1). This is in sharp contrast with the nontrivial 
geometric phase factor of —1 for the case of intermedi- 
ate electron correlations (2 < ?7 $ 4) . A reentrant behav- 
ior of the trivial geometric phase factor of -1-1 occurs for 
the case of strong antiferromagnetic electron correlations 



We examine the exchange symmetry of polarons in dif- 
ferent aspect. For the case of weak electron correlations 
(or free electron limit, U ^1) the predicted exchange sym- 
metry indicates a bosonic nature of polarons with the ab- 
sence of fictitious magnetic flux line. This is depicted in 
Fig. @(a) . As correlation (repulsive interaction) between 
electrons increases, a change of the geometric phase is 
predicted. That is, for the case of intermediate correla- 
tions it can be stated that the polarons hop around a ficti- 
tious flux of an elementary flux unit (fluxon) , thereby ex- 
hibiting a fermionic nature (see a solid line in Fig. p(b)). 
For the case of strong antiferromagnetic electron corre- 
lations, the exchange of the two polarons indicates the 
presence of two fictitious flux quanta, as is depicted in 
Fig. ra(c). Then the statistical transmutation of boson 
to a composite boson occurs. By composite boson we 
mean that the boson remains as a boson with the ac- 
quired phase of 2tt. In short the number of fictitious 
flux quanta is seen to vary with the strength of electron 
correlation and allows the statistical transmutation. 






FIG. 6. Exchange of two polarons. Open circles denote 
polarons, and each vertical solid line represents a fictitious 
flux line of one flux quantum (flux unit). 



IV. CONCLUSION 

By applying the Lanczos exact diagonalization method 
to Holstein-Hubbard, Holstein-tJ, and Holstcin-tJ^ 
model Hamiltonians, we have examined the different as- 
pects of geometric phases acquired by electronic wave 
functions during the transport of a single polaron around 
a triangular loop. We find that the nontrivial geometric 
phase factor of —1 is caused by the presence of trans- 
verse spin fluctuations. Further by examining the ex- 
change symmetry of polarons based on the Holstein- 
Hubbard model we discover the statistical transmutation 
of polarons which varies with the strength of antiferro- 
magnetic electron correlation. It is interesting to note 
that polarons behave as bosons in the hole-doped two- 
dimensional system (background) of weakly and strongly 
correlated electrons and as fermions in the background 
of intermediate correlation strength. In the future it re- 
mains to examine how the statistical transmutation for 
the exchange of polarons occurs as a function of electron 
correlation and whether there exists a possibility of any- 
onic exchange phase other than 6 — and 6 — tt, that is, 
< S < n. 
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APPENDIX: 

For the case of multi- valued real wave functions here we 
prove the Berry phase to be 7(C) = 6'(R(T))-6'(R(0)) by 
introducing a local gauge transformation in the parame- 
ter space, R(= u). The Berry phase 7„(C) is a geometric 
contribution to the total phase change for the final state 

of mr)) = eMnniC)]e^p[-j- J^ dtE^iRit))} j^-W) 



where T is the time period. El For a subsystem of fast de- 
grees of freedom coupled to a subsystem of slow degrees 
of freedom the Berry phase is obtained from 



7„(C) =1 /(n;R|Vn|7i;R)-dR 



(Al) 



as is well known. Q 

The equation of motion for the subsystem of fast de- 
grees of freedom is given by 



H(R{t))\n; R(i)> - ^„(R(i))|n; R(i)) 
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where R(i) is the slowly varying parameter. To obtain 
the nontrivial (non-zero) Berry phase in Eq. (Al), the 
eigenstate |7i;R(i)) should be complex, non-degenerate, 
and single- valued. The gauge potential 



A(R) = i{n;'R\Va\n;R) 



(A3) 



vanishes for the case of the real eigenstates |n; R). 

It is obvious that the nontrivial (nonzero) 'magnetic 
field' B(R) = V X A(R) cannot be defined for the case 
of real eigenstates \n;'R,{t))r, since the gauge potential 
A(R) vanishes. However, by using the following gauge 
transformation it is possible to define the nonvanishing 
Berry phase: 



;R) 



= e-'''(^)|n;R) 
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Indeed, this is the local gauge transformation in the pa- 
rameter space R, which provide s th e gauge invariance of 
the Schrodinger equation, Eq. (A2). Using Eq. (A4) in 
Eq. (|A1| ) we readily obtain the Berry phase 7„(C), 



7„(C) = 0(R(r))-0(R(O)) 



(A5) 



or usingthc notation used for the parameter space in 
Section y, 

-/n{c) = eiuiT)) - eiuio)) . 
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